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Abstract: We estimate the bounds on Higgs and sparticle masses and discuss their cor- 
relations in the constrained minimal super symmetric standard model (CMSSM). In our 
analysis we have applied the present constraints from collider and low energy experiments, 
as well as the experimental bound on cold dark matter from WMAP. For a given lightest 
Higgs boson mass, which is expected to be measured with good precision at the LHC, we 
find important correlations between the Higgs and sparticle masses which allows one to 
delineate the MSSM model parameters and particle spectra. We have also demonstrated 
an important complementarity between the LHC and direct dark matter detection experi- 
ments emphasizing that by including the experimental input both from collider physics and 
from dark matter detection experiments, one would significantly improve the measurement 
of the SUSY spectrum and the underlying parameter space. 
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1. Introduction 



In recent years the measured value of the top quark mass has gradually inched its way 
down to its current value of 170.9 db 1.8 GeV pj, which is significantly below the central 
value of around 178 GeV used in various studies just a few years ago. This downward shift 
of the top quark mass has important implications for the particle spectrum of the minimal 
supersymmetric standard model (MSSM), in particular for the lightest CP-even Higgs 
boson mass, rrih, since the leading radiative corrections to rrih are proportional to the square 
of the top mass, mf , while rrih > 114.4 GeV from LEP2 in the decoupling limit @. At tree 
level rrih < Mz\ cos 2(3 1, and so significant radiative corrections are necessary to overcome 
the LEP2 bound. In the absence of large trilinear couplings, significant radiative corrections 
can be generated only if the stop masses are in the TeV range. Alternatively, significant 
stop mixing in the presence of large trilinear couplings can drive rrih up without excessively 
heavy stops. In this paper, we investigate both possibilities and study the correlations of 
the Higgs and sparticle mass spectrum in the constrained MSSM (CMSSM) ^ [|, [|, |(| |?J 
under the plausible assumption that the neutralino constitutes the cold dark matter in the 
universe. 

The CMSSM has received a great deal of attention in recent years, especially for m% 
close to 178 GeV J§ |, 0, |Tl|, 0, [13|, [[§ [TjJ g]| [H], 0. When we started this analysis, 
the central value was writ = 171.4 ± 2.1 GeV [IS], which has recently drifted slightly lower 
to 170.9 ± 1.8 GeV Q. Our results are based on rrit = 171.4 GeV. Note that we have found 
agreement between our results with partially overlapping results of some recent papers 
which also use an updated value for the top-quark mass 21, ^2[ |23[| . 

The CMSSM contains the following five fundamental parameters (we follow notations 
and conventions of Ref . pl| ) 



m , m 1 / 2 , A , tan/3, sign(/x) , 



where wiq is the universal soft supersymmetry breaking (SSB) scalar mass, rrii/2 the uni- 
versal SSB gaugino mass, and Aq the universal SSB trilinear scalar interaction (with the 
corresponding Yukawa coupling factored out). The values for these three parameters are 
prescribed at the GUT scale, Mqut — 2 x 10 16 GeV. tan (3 is the ratio of the vacuum 
expectation values (VEV) of the two Higgs doublets at the weak scale, and /i is the super- 
symmetric bilinear Higgs parameter whose magnitude, but not the sign, is determined by 
the radiative electroweak breaking conditions 1 . 

The results for CMSSM studies are usually presented in the (mo,mi/ 2 ) plane. How- 
ever, anticipating that the lightest MSSM Higgs will be found at the LHC and that its 



mass according to ATLAS [28| and CMS 2S] technical design reports, will be determined 
to a precision which is better then 1%, it is worthwhile to understand possible correlations 



historically this framework is also called mSUGRA and is associated with supergravity. However, 
supergravity does not necessarily lead to high scale universality as originally thought. Occasionally, the 
abbreviation mSUGRA is used only for a subset of ( [l.l| ) where some additional relations coming from a 
particular form of Kahler potential are assumed (see e.g. p6 ). To avoid possible confusion, we decided 
to use the term CMSSM and refer interested readers to Ref. [£4, ^tJ for textbook discussions. 
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between such a precisely measured Higgs mass and sparticles masses as well as the un- 
derlying model parameters. Such correlations would allow one to delineate the underlying 
SUSY model and, moreover, to predict some model parameters and masses, such as the 
trilinear coupling, At, whose measurement at the LHC may be problematic. The mass rrih 
may be experimentally determined to an accuracy better than a GeV, which is better than 
the theoretical calculation of at present J3C[ |. In this study, we assumed a theoretical 
uncertainty Arrih of 1 GeV in estimating the precision with which the sparticles masses 
and model parameters are predicted as a function of m^, keeping in mind the optimistic 
hope that the theoretical uncertainty for the Higgs mass calculation will be improved by 
the time of the actual Higgs boson search at the LHC. In our analysis we require that the 
magnitudes of m^j and Aq should not exceed 2 TeV, while mo is taken to be < 5 TeV. 
This choice covers the whole dark matter-motivated CMSSM parameter space except the 
focus point region (see discussion below), which can go with tuq all the way up to Mqut, 
and turning the CMSSM scenario into a split-SUSY one [31]. 

The sign of the MSSM parameter fi is taken to be positive, as dictated by the experi- 
mental measurements of the anomalous magnetic moment of the muon, Aa^. We take into 
account the LEP2 constraints on the chargino mass, experimental data on rare B decays, 
and the WMAP relic density constraints. 

The paper is organized as follows. In Section |2| we review details of the MSSM pa- 
rameter space scan and implementation of the various constraints. In Section |3] we present 
our results on various correlations of the MSSM spectra which would be crucial for LHC 
collider phenomenology, and in Section ^| we summarize our conclusions. 



2. SUSY constraints and scanning procedure 

For the MSSM spectrum calculation we use the IsaSUGRA program, which is a part of 
IS A JET 7.74 package [p2]| . In this package, the weak scale values of gauge and third 
generation Yukawa couplings are evolved to Mqut via the MSSM renormalization group 
equations (RGEs) in the DR regularization scheme, where Mqut is defined to be the scale 
at which g\ = 52 • We do not enforce an exact unification of the strong coupling 53 = 51 = 52 
at Mqut, since a few percent deviation from unification can be assigned to unknown GUT- 



scale threshold corrections [33, 34]. At Mqut, the universal soft supersymmetry breaking 
(SSB) boundary conditions are imposed and all SSB parameters along with the gauge and 
Yukawa couplings are evolved back to the weak scale Mz- In the evaluation of Yukawa 
couplings the SUSY threshold corrections [35] are taken into account at the common scale 
Msusy = \J m t L m t R > an d the entire parameter set is iteratively run between Mz and Mqut 
using full 2-loop RGEs until a stable solution is obtained. To better account for leading-log 
corrections, one-loop step-beta functions are adopted for gauge and Yukawa couplings and 
the SSB parameters mt are extracted from RGEs at multiple scales m,j = rrii (mi). 2 The 
RGE- improved 1-loop effective potential is minimized at an optimized scale Msjjsy, which 
effectively accounts for the leading 2-loop corrections. Full 1-loop radiative corrections are 



2 Detailed discussion of SUSY threshold effects in IS A JET and other public codes can be found in ref. [B6| 
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incorporated for all sparticle masses according to ref. [35]. The numerical results are in 
close accord with those generated by the FeynHiggsFast program 



The requirement of radiative electroweak symmetry breaking (REWSB) |38|, 40 
|4l|| puts an important theoretical constraint on the parameter space. Another important 
constraint comes from limits on the cosmological abundance of stable charged particles [42] 
- this excludes the parameter points where charged SUSY particles, such as fi or t\, become 
the lightest supersymmetric particle (LSP). 

We also impose bounds on sparticle masses resulting from direct searches at colliders. 
The LEP2 search for pair production of sleptons [43] puts the lower limit mr > 99 GeV 
provided the mass difference mj — > 10 GeV. Another important LEP2 constraint 
comes from the lightest chargino Q|, th^ > 103.5 GeV. 



Recent search results for squark/gluino production at the Tevatron [Ha, 46] put bounds 



on the gluino mass m g > 370 GeV for rriq ~ rrig, and m g > 220 GeV for rriq < rrig. In 
models with gaugino mass universality, like CMSSM, m g ~ 3.5m^. Thus, the LEP2 
chargino bound translates into a bound on the gluino mass, m g > 350 GeV, to which 
Tevatron comes close only for rriq ~ m g . 



Tevatron searches for direct pair production of stops |t7], [48| impose the bound > 



132 GeV for < 48 GeV. However, the requirements of GUT-scale gaugino mass uni- 
versality tell us that ~ 0.5m^ > 51 GeV for points satisfying the LEP2 chargino 
bound. This chargino- neutralino mass ratio gets smaller in the HB/FP region at large mo, 
but here the stops are always heavier than about 500 GeV, as we will show later. 

Similar searches for direct pair production by DO collaboration [49| put a bound > 
222 GeV for m% < 90 GeV; the CDF limit is slightly lower, rn~ bi > 193 GeV for < 
80 GeV p7[]. In CMSSM, the lightest sbottom is dominantly left-handed whose mass is 
given by the approximate relation m~ ~ 0.52m.Q + 5m^ 2 . The lightest chargino is mostly 
wino and obeys the approximate relation m^ ~ |"ii/2- Therefore, the LEP2 chargino 
bound leads to m^ > 350 GeV, significantly above the Tevatron limit. 

The LEP2 Higgs bound || further constrains the MSSM parameter space. After the 
LEP2 direct SUSY searches constraints are applied, the pseudoscalar mass rriA is limited to 
be above 200 GeV (as we present later on), which ensures the decoupling regime and SM- 
like coupling of h to gauge bosons. Therefore, the SM-Higgs LEP2 bound rrih > 114.4 GeV 
is applicable to the CMSSM scenario. This bound can be very different in the MSSM light 
Higgs boson scenario when the LEP2 bound allows a higgs scalar as light as about 
60 GeV. We would like to point out that in our plots we do not remove points with the 
Higgs boson mass below the SM-higgs LEP2 bound, but indicate its position whenever 
possible for the sake of clarity. 

A combination of the recent WMAP data with other cosmological observations leads to 
tight constraints on the cold dark matter (CDM) relic abundance, and the most important 
and dramatic constraint on the CMSSM parameter space. The exact numerical value 
depends on a number of assumptions about the history of the early Universe 51] as well as 
on a combination of data sets chosen. In our analysis we use the constraints arising from 
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a combination of WMAP and the Sloan Digital Sky Survey data [g^] 

^cdm/i 2 = 0.111±g;g^ (2a) , (2.1) 

which assumes standard ACDM cosmology. We further assume that the bulk of CDM 
is composed of the lightest neutralino Z\ which was in thermal equilibrium in the early 
Universe 3 . To evaluate the neutralino relic density we employed the IsaReD code |54| (part 
of IsaTools package) , which uses several thousands 2 — ► 2 annihilation and co-annihilation 
Feynman diagrams generated by the CompHEP package pq , [56| . 

An important constraint comes from the measurement of the muon anomalous mag- 
netic moment a M = ^ ^ by the E821 collaboration |57|| . The current experimental value, 
af p = (11659208.0 ± 6.3) x 10~ 10 deviates by 3Aa from the SM prediction which uses 
hadronic vacuum polarization determined from e + e~ annihilation data, but there is no sig- 
nificant deviation from the SM predictions if the hadronic contribution is calculated using 



r-decay data |58|, |59|, |60|]. There is a stronger and stronger tendency at present to prefer 
e + e~ data to evaluate the lowest-order hadronic contribution to a„. The most recent analy- 
sis performed in ref. |p8| , based on e + e~ data, yields a deviation Aa^ = (27.5 ±8.4) x 10~ 10 
from the experimental value. In view of the lack of consensus on the computation of the 
SM part [|58], |5^, ^0|, we have decided not to apply the (g — 2) M constraint but to indicate 
the 3a allowed region for (g — 2) M , 

3.4 x 10~ 10 < Aa M < 55.6 x 10~ 10 , (2.2) 

with special colors to allow readers to make their own judgment on this constraint. This 
3a interval uses (g — 2)^ predictions based on e + e~ data and was taken from the latest 
summary of (g — 2) M status |h| . 

We used IsaAMU subroutine described in ref. to evaluate the SUSY contribution 
to (g — 2)^. The contribution of SUSY particles comes from WiU^ and Zifli^ loops, where 
the total chargino contributions dominates in CMSSM. Therefore, (g — 2) M puts an upper 
limit on the masses of charginos, neutralinos and second generation sleptons. For tan (3 
not too small, Aa^ USY cx [iM.2 tan/3/m 4 , where M2 is the SU(2) gaugino mass and fh is 
the heaviest mass in the loop. Thus, it is particularly important for large tan/3 values and 
light sleptons and -inos (i.e., small mo and mi/2 region). Note also that for models with a 
positive gaugino mass, like CMSSM, the sign of \x has to be positive. 

The good agreement of the experimental value of the flavor changing b — > 57 decay 
rate with the SM result imposes an important constraint on any beyond the SM physics. 
In MSSM, the most important SUSY contributions come from tH^ and t loops. To 



evaluate the SM+MSSM branching fraction, we employed the IsaBSG subroutine |6^, p4| , 
which uses the effective theory approach. Namely, when any sparticle threshold is crossed, 
the sparticle is integrated out, thereby inducing a new effective-operator basis. The re- 
sulting Wilson coefficients are evolved down to Q = mt,, thus summing large logarithms 
that can occur from the disparity between the scales at which various particles enter the 
loop corrections. In the CMSSM, the assumed universality of soft SUSY breaking (SSB) 



s Non-thermalized neutralinos were recently discussed in ref. |53 
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masses leads to a particularly simple framework known as minimal flavor violation (MFV) 
in which flavor mixing arises only from the CKM matrix. The experimental world average 



measurement for Br{b — > 57) is known to be (3.55 ±0.26) x 10~ 4 [65], while its updated SM 
prediction has been recently shifted down to (3.15±0.23) x 10~ 4 value [66, 67]. Combining 
the experimental and theoretical errors in quadratures we apply the following constraints 
at 2cj level in our study: 

2.85 x 1(T 4 < Br(b -» S7) < 4.24 x 10~ 4 . (2.3) 

This conservative approach has two advantages. Firstly, it allows one to accommodate 
the theoretical uncertainties that mainly come from the residual scale dependence and are 



~ 12% [ 38 1 . Secondly, the SSB terms might have small off-diagonal entries at Mqut that 
will change Br(b — > 57) without significantly affecting other quantities. 

Another B-physics constraint comes from the upper limit on B s — > ^ + /^ _ decay branch- 
ing fraction, 

BF(B S -> < 1.0 x 10~ 7 , (2.4) 



as was reported by the CDF collaboration 6£j. We compute the SUSY contribution us- 



ing IsaBMU code jTOj ] which assumes MFV. The SM predicts a very small value for this 
branching fraction, namely BFsm(B s — > ju + /i~) ~ 3.4 x 10~ 9 , while the SUSY contribution 
behaves as tan 6 f3/m\ and hence is particularly important at large tan j3. However, it turns 
out that this constraint is always superseded in the CMSSM parameter space under study 
by the one from b — > S7. 

To find bounds and correlations for the Higgs and sparticle masses, we performed 
random scans for the following range of the CMSSM parameters: 

< m < 5 TeV, < m 1/2 < 2TeV 

A) = 0.5TeV, 0, -ITeV, -2TeV, tan/3 = 5, 10, 50 and 53 (2.5) 

with jj, > and mt = 171.4 GeV. We have also performed scans for intermediate values 
of tan j3 = 30, 45, but found results to be similar to tan (3 = 10, 50 cases respectively 
and decided to omit them. The use of IsaTools package for implementation of the various 
constraints mentioned above was crucial in our study. 



3. Results 



Let us first demonstrate our procedure by consecutively applying the constraints mentioned 
in the previous section to the case with Aq = —2 TeV and tan/3 = 10. We start by 
presenting our results in the conventional (m-o, ^1/2) plane illustrated in Fig. |]. Frame a) 
displays the colored points allowed by radiative electroweak symmetry breaking (REWSB) 
and from the requirement that the neutralino is the LSP 4 . 

4 In fact, to keep the data set size manageable, we show only points where Q.h 2 < 10. The CMSSM 
points with too high a relic density, denoted by gray or white color above the gray-colored area, are not 
relevant to the current discussion. 
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tan$=10, A =-2TeV,\i>0 




il/T< 0.122 
b-^sy excl. 



\ 




(g-2) u within 3cs 
(g-2)^ outside 3a 




2 4 2 4 2 4 

m (TeV) m (TeV) m (TeV) 

Figure 1: (mo,m 1 / 2 ) plane for tan/3 = 10, Aq = —2 TeV. In frame a), the gray regions are 
excluded by WMAP, blue ones are excluded by LEP2 chargino bound, and the yellow region has 
m,h below 114.4 GeV; the remaining colors denote various WMAP-allowed regions described in the 
text. In frame b), the green regions satisfy both the WMAP and b — > sj constraints, while the red 
region satisfies the WMAP bound but fails the b — > constraint. In frame c), regions surviving 
the WMAP and b — ► 57 constraint are presented; light (dark) green region corresponds to (g — 2) M 



within (outside) the 3cr range (2.2) 



The gray region at low mo and large-to-medium values at the left is excluded be- 
cause m-fj < nig , while the white bulge at low mo and low myj under the t-coannihilation 
contour is excluded because the top-squark is LSP in this region (typical for such a large 
negative Aq). Finally, the white corner region at the bottom-right below the colored band 
is excluded because of the failure of the REWSB condition. The LEP2 constraint on 
chargino mass removes the points shown in blue color. After application of the WMAP 
upper bound on LSP DM, the allowed region is drastically shrunk to several small regions 
of the parameter space, with the remaining gray points associated with unacceptably high 
DM relic density. This exhibits a well known fact that in CMSSM the neutralinos do not 
efficiently annihilate in the early universe except in a few very special narrow regions of 
the parameter space 5 : 

• the bulk region at low tuq and low TO1/2, shown in brown color, where neutralinos 
efficiently annihilate via t— channel slepton exchange. This area is almost completely 



ruled out by LEP2 searches 0> El L3> Ell 



These regions sometimes appear discontinuous on some plots only due to the finite resolution of our 
scan. 
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the stau co-annihilation region (fZ) j75], 76, 77, [78|, [7|, 8C, |82], 54] at low mo, 
shown in magenta color, where — m^ and rapid reactions Zifi — » X and 

fxri — > X (here X stands for any allowed final state of SM and/or higgs particles) in 
the early universe lower the neutralino relic density to the WMAP range (|2.1|); 



the stop coannihilation region (iZ) appearing at low thq and low m^2, but for par- 
ticular values of A$, shown in orange color, where rug — [83, |4|, 85]; 

the hyperbolic branch/focus point region (HB/FP) at large mo, shown in turquoise, 
where \fx\ is small and the neutralino develops a substantial higgsino component that 



enhances its annihilation into WW, ZZ and Zh pairs in the early universe [86, 87 



, 89, 90|; 



the h— corridor (HF) at low mi/2; shown in purple color, where 2m % — and 
efficient annihilation through the light MSSM Higgs boson resonance occurs in the 
early universe |9l], [73| ; 

the A— annihilation funnel (AF) that appears only at suitably large values of tan (3 
and is therefore absent in the current plot. In this case, 2mj; ~ m^ and resonance 
annihilation of neutralinos via the broad A and H Higgs bosons becomes possible 



93j, |94|,|95|, |96|,|9J. 



Next, in frame b), we show the effect of the b —* sj constraint on the WMAP allowed 
region: red points denote the excluded region where the branching fraction exceeds the 



range (2.3) which covers the bulk region, stop coannihilation and the low m^ part of the 
stau coannihilation regions. The green points survive the combination of the WMAP and 
b — > S7 and are the subject of our analysis in this paper. Finally, in frame c) we illustrate 
the effect of — constraint fl2.2j ) which, at 3a level, excludes the h-corridor and HB/FP 
regions at medium and high mo, denoted by dark green color. This occurs since sleptons are 
heavy in these regions and their loop contributions do not make significant contributions 
to {g — 2)^. On the contrary, in the fZ region, shown in light green color, the sleptons are 
light and can account for the observed deviation ( |2.2| ). 

In Fig. ^ we present results for the same parameter space, but in the (m^m/i) plane 
with the color coding as in Fig. [l|. Note that the LEP2 Higgs boson limit excludes the tZ re- 
gion, leaving fZ, HF and HB/FP regions respecting the WMAP dark matter constraint. 
It is very noticeable, that the regions allowed by SUSY constraints appear as very narrow 
bands exhibiting an obvious correlation of tjia with m^. This correlation is related to the 
two- fold solutions mainly defined by the DM constraints: the lower band corresponding to 
fZ region (connected to tZ and AF regions at high values of tan/3) and the upper band 
corresponding to the HB/FP region (and connected to HF region at low to intermediate 
values of tan (3) . 

Qualitatively similar patterns can be observed for other Aq and tan (3 values and are 
presented in Fig. |3[ Here gray, light green, light blue, and orange colors correspond to 
tan/3 =5, 10, 50 and 53 respectively and satisfy the 3a bound on Aa^ given by Eq. . 
Black, dark green, blue and red colors correspond to the region outside the Aa M 3a bound; 



-7- 



tan$=10, A =-2TeV,\i>0 




110 120 110 120 110 120 

m h (GeV) m h (GeV) m h (GeV) 

Figure 2: (mf,, m^) plane for tan (5 — 10, = —2 TeV. The color code is the same as in Fig. |l|. 



frames a),b,)c),d) present results for Aq = 0.5, 0,-1 and —2 TeV respectively. The upper 
and lower bands for each color represent the HB/FP(+HF) and fZ{+tZ+AF) regions 
respectively. One finds that the CP-odd Higgs boson mass is always above ~ 200 GeV, 
corresponding to the decoupling regime and SM-like nature of h. One can see that for 
Aq = 0.5 TeV and (frames a) and b) respectively) the parameter space for tan/3 = 5 
is completely excluded by the LEP2 limit, > 114.4 GeV. This is due to the low tree- 
level Higgs mass and the absence of large enough radiative corrections. However, for 
Aq = — 1 TeV and —2 TeV (frames c) and d) respectively), enhanced radiative corrections 
to the Higgs boson mass open up some viable CMSSM parameter space even for tan (3 as 
low as 5. We note here that at the EW scale mh depends, to a good approximation fl98| , 
on the absolute value of Aq. Starting with negative values of At at MquTi it turns out 
that At evolves to even larger negative values at the EW scale, thereby enhancing mh and 
relaxing the lower bounds on the masses of the SUSY particles. This explains our choice 
for Aq values used in the paper. 

We stress that the light colors in Fig. || and subsequent figures represent the parameter 
space which respect all of the constraints mentioned above including {g — 2) M , while the 
dark colors represent the parameter space respecting all constraints except (g — 2)^. For low 
and intermediate values of tan (3 only some portion of the f-coannihilation regions (lower 
curves) denoted by gray (tan (3 = 5) and light green (tan (3 = 10) colors satisfy all of the 
constraints, while the focus point region is excluded. In fact, this known feature of CMSSM 
represents the tension between essentially two constraints: b — > S7 and {g — 2)^. The root 
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Figure 3: Allowed region for CP-odd Higgs boson mass tua versus mh- Gray, light green, light 
blue, and orange correspond to tan f3 =5, 10, 50 and 53 respectively and satisfy the 3a bound on 



Aa^ (2.2). Black, dark green, blue and red correspond to tan/? =5, 10, 50 and 53, with Aa p outside 
the 3cr range. 



cause of this tension is the universality condition in CMSSM which implies equal squark 
and slepton soft mass parameters at Mqtjt- On the other hand, b — > 57 data agrees with 
the SM prediction (implying heavy third generation of squarks at the EW scale) whereas 
there seems to be quite a significant deviation between the SM prediction and (g — 2)^ 
measurement (implying light second generation of sleptons at the EW scale). It has been 
shown [15] that this tension can be resolved either with large tan /3 values or, alternatively, 
by relaxing the GUT-scale generation universality condition. Indeed, one can see that at 
large tan (3 there are allowed regions with both upper and lower (light blue (tan (3 = 50) 
and orange (tan/3 = 53)) curves. In the case of large tan/3, tua (run) up to about 1 TeV 
can be accessible at the LHC and used for a consistency check of the CMSSM. 

In Fig. |I| we present correlations between the gluino and the light Higgs boson masses 
with the same color coding as in the previous figure. The pattern of bands in the (mg, m^) 
plane is quite different from the one for (mA,mh) in the previous figure. Gluino mass is 
mainly defined by m^ 2 which, in turn, is similarly correlated with along the fZ and 
HB/FP regions. In the (mg, m^) plane the bands are close to each other. The gluino 
mass increases in the 'far' HB/FP region of large mo and large rn^j, since in this region 
m i/2 goes beyond 1 TeV where the fZ region typically ends. The 'far' HB/FP region is 
characterized by the heaviest and heaviest gluino mass in this particular plane. Notice 
that the gluino mass is always above 400-500 GeV, so that the recent Tevatron bounds on 
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gluino mass [45, 46] do not affect the parameter space left after imposing DM, LEP2 and 
b — ► S7 constraints. On the other hand, the LHC with 100 fb _1 luminosity will be able to 
probe the parameter space corresponding to m~ g up to ~ 3 TeV [88. [9(], ^9|, |100| , |101| , |102| , 
covering almost the entire fZ band. Also, the LHC with 100 fb _1 can indirectly probe rrig 
up to ~ 2.5 TeV in the HB/FP region via observation of lepton + jets+ $t signal from 
lighter gauginos [|103|| . Distinguishing fZ and HB/FP bands in (m^m^) plane could be 
somewhat problematic, even though, potentially, the knowledge of the light Higgs boson 
mass at the percentage level would allow one to indirectly determine the mass of very heavy 
gluinos with 5 — 10% precision. 

In Figs. H and |6] we present plots for tri^ versus and versus rrih, respec- 
tively. The upper band of the allowed parameter space for each particular value of 
tan/3 corresponds to the fZ(-\-tZ+AF) region. In this region, the mass of the chargino 
and the neutralino is being driven up with the rise of tan/3 up to about m^y — 1.7 TeV 
and ~ 850 GeV for tan/3 = 53. The HB/FP and HF regions are characterized by 

small |/z| values resulting in light quasi-degenerate Z\ and W\ states. This happens as 
a result of the following non-trivial interplay of the RG equations. With m^ 2 fixed, in- 
creasing mo leads to the suppression of terms with top and bottom Yukawa couplings in 
the RG equations for rriH u and rriH d - This effect drives the higgs SSB masses towards the 
"no REWSB" condition and respectively lowers the value of [//[. These regions form the 
lower band of the allowed parameter space. The upper tip of this band can reach the mass 
region as large as ~ 700 GeV. This limit is defined by the maximal value of rriQ = 5 TeV 
in our scan and practically does not depend on tan (3. One can see that the slope of the 
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Figure 5: Allowed region for the lightest chargino mass versus m/,. The color code is the same 
as in Fig. || Dashed (dotted) horizontal lines represent approximate reach of ILC500 (ILC1000). 
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coannihilation band is typically larger than the slope of the allowed band of the parameter 
space corresponding to the focus point region. Therefore, in the coannihilation region one 
can expect better accuracy in the chargino/neutralino mass determination through the 
correlations presented in Figs. [| and |(| Also, using rrig — rrih correlation, when the gluino 
is observed and reconstructed, one can actually predict the chargino and the neutralino 
masses, whose observation in EW production processes could be problematic. To put things 
in perspective, we show in Fig. [| the kinematic limits for W\W\ pair production, which 
can serve as good approximations for the reach of ILC500 (yfs = 500 GeV) and ILC1000 
(y/s = 1000 GeV) machines [104, 105, 10C]. We see that the linear collider can access 
only the lower end of the coannihilation band (corresponding to low m\/2 values), while 
the HB/FP region can be probed almost entirely. This is in contrast to the LHC, which 
has large coverage in the coannihilation region, but can only probe the lower half of the 
HB/FP region, as was shown in Fig. 0. 
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Figure 7: Allowed region for top squark mass mass versus to/j. The color code is the same as in 
Fig. |. 

In Fig. [j] we present correlations between the lightest top-squark mass and m^. In this 
figure, the upper curves for each value of tan/3 correspond to the HB/FP(+HF) region 
and the lower curve to the ?Z(+tZ+KF) region. At tan/3 = 53 and Aq = 0.5 TeV and 
Aq = in frames a) and b) respectively, one can see that the orange curve is widened 
in its upper part. This region corresponds to the AF region, where the large range of 
radiative corrections to the top squark mass and its stronger dependence on various model 
parameters spoil somewhat the — correlation. In this region, mi > 500 GeV, 
thus exceeding the Tevatron bound 47, 45]. One can also notice that for Aq = — 1 TeV 
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and Aq = — 2 TeV, presented in frames c) and d), respectively, the top squark exhibiting 
stop-coannihilation region is allowed to be as light as 200 GeV. This region is represented 
by the lower tip of the light-green line in frame c), and the lower tips of gray and light- 
green line in frame d). As Aq becomes increasingly negative, the mixing between and 
tR increases (thereby reducing the mass of t\) which, in turn, drives to larger values 
through radiative corrections. This explains why the curves shift to the lower right corner, 
as we go from frame a) to d). 

Let's take a closer look at the 'anatomy' of the stop co-annihilation region. In Fig. || 
we present constraints similar to Figs. |l] and ^, but in the (m^, m/j) plane with the same 
color coding for tan/3 = 30 and Aq = — 1 TeV. The DM-allowed area consists of HB/FP 
region (upper branch), fZ (lower branch), and tZ (the very bottom tip of the DM allowed 
region) regions. From the frame b) one can clearly notice that the most serious constraint 
for the stop-coannihilation region comes from the b — > s'y. Indeed, for a light top squark, 
the contribution from the stop-chargino loop increases Br{b — > 57) beyond the acceptable 
level. The reason for the appearance of the very small allowed island surrounded by the 
b — > 57 excluded region is the non-trivial cancellations between stop-chargino and top- 
charged Higgs loops. Thus, the tZ region does survive in the CMSSM framework, but this 
scenario is highly constrained. 
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Figure 8: Allowed region for top squark mass mass versus to/j. The color code is the same as in 
Fig. 0. 

Next, in Figs. |9| and we present correlations between the bottom squark mass and 
mh, and stau mass and mh, respectively. One can observe qualitatively similar pattern of 
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correlations in these two figures, where again, as in Fig. [7|, the upper curve corresponds 
to the HB/FP region, while the lower curve corresponds to the stau-coannihilation one. 
At large tan (5 and Aq = 0.5 TeV and Aq = 0, the A— funnel region opens up, which is 
reflected in the behavior of the orange curve. Increasing tan (3 boosts up the bottom and 
tau Yukawa couplings, which reduces the sbottom and stau SSB masses through RGE 
effects, and increases L-R mixing; both effects tend to reduce b\ and f\ masses. This effect 
is more noticeable for staus, as seen from the wider separation of the HB/FP bands in 
Fig. [l(]. Notice that in Fig. [9|, the minimal value of bottom squark mass is ~ 500 GeV, 
which is significantly above the current Tevatron bound [fl?], |2^| . Kinematical limit for 
W\W\ pair production at ^fs = 1000 GeV, shown by the dashed horizontal line, indicates 



an approximate reach of ILC1000 [|10j, [L0|, [UKf . 

It is quite informative to also consider the correlation plots for a fixed value of tan (3 
and different Aq values in the same frame. For example, in Fig. 11 we present the allowed 
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Figure 11: Allowed region for CP-odd Higgs boson mass versus m^. Gray, light green, light blue, 
and orange correspond to Aq = 0.5, 0,-1 and —2 TeV, respectively and satisfy 3er bound on Aa M 



(2.2). Black, dark green, blue and red correspond to A = 0.5, 0, —1 and —2 TeV and have Aa M 
outside the 3a range. 



region for the CP-odd Higgs mass versus m^. Gray, light green, light blue, and orange 
regions correspond to Aq = 0.5, 0, —1 and —2 TeV respectively and satisfy the 3<r bound 
(|2,2| ) on Aa M ; black, dark green, blue and red colors correspond to regions outside the 
Aa^ bound. If we assume that tan/3 is measured, which can be done at the LHC with 
an accuracy better than 20% for low to intermediate values of txia |28[| , then the strong 
rrtA — correlations above can be used for an indirect prediction of uia, in case the 
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CP-odd Higgs boson is too heavy to be directly accessible. A 1% accuracy in the 
measurement would allow one to estimate uia with a precision of 5-10%. 

One can also use these correlations to determine (or at least strongly constrain) the 
trilinear coupling Aq, whose direct measurements are quite problematic. In Fig. [12], we 
present correlations between the gluino and light higgs boson masses, where we group vari- 
ous Aq values for fixed tan (5 in each frame. For each Aq value, the upper band corresponds 
to the coannihilation region and the lower one to the HB/FP region. The two groups of 
bands are very close, but if we discover that we are in the HB/FP or in the coannihilation 
region, half of the curves will be removed and the remaining ones are well separated for 
tan/3 > 10, which can be used to extract the Aq value. 
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Figure 12: Allowed region for gluino mass versus m/,. The color code is the same as in Fig. |Tl| . 

Furthermore, the study of the correlations between various SUSY masses would allow 
one to further delineate the SUSY parameter space. For example, in Fig. [l^ we plot the 
allowed region in (m^ , m^) plane. We see that the two-fold solution bands form two well 
separated groups - the lower one representing the coannihilation region and the upper one 
corresponding to the HB/FP region. One should also notice a third type of bands, namely 
the vertical one, located at low values and corresponding to the HF region. Given 
rriA (either from direct measurements or extracted from Fig. 11) and tan/3 values, one 
can deduce Aq from the reconstructed neutralino mass. One should point out, however, 
that the LHC potential can be quite limited in the reconstruction of the heavy neutralino 
masses in certain regions of the parameter space |102[| . Analyzing the direct DM detection 
(DD) rates, we have found an important complementarity of DM detection experiments in 
this respect, since these experiments can cover a much wider range of neutralino masses. 
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In Fig. [l4|, we present the DD rates, in detector-invariant way, as spin-independent 
neutralino-proton elastic scattering cross section versus the neutralino mass. We used 
IsaReS code | 107 |, a part of IsaTools package, with pion-nucleon £ term assumed to be 
45 MeV 6 . One can clearly read patterns coming from different dark matter motivated 
regions. The upper band corresponds to HB/FP region, where Z\ has a substantial higgsino 
component thereby enhancing its scattering cross section off the proton. The lower pattern 
comes from the stau coannihilation region (merging with A- funnel for tan/3 = 53). We 
see that the HB/FP and coannihilation regions are always separated in terms of spin- 
independent cross section. One should also note the h-funnel region, which is represented 
by narrow vertical band at ~ 55 GeV. This region is well separated from HB/FP and 
fZ regions and could be almost completely covered by XENON 1 ton detector. 

The current best limit comes from the XENON-10 collaboration [10S], which obtained 
an upper limit a si {Zip) £ 8 x 1CT 8 pb for trig ~ 100 GeV. We also show the projec- 
tion for CDMS2 |110| ], which is expected to release results next year, and for its planned 
upgrade, SuperCDMS [111]. The reach of XENON-1 ton detector as representative of 



many planned large noble gas detectors [112, |113] , 114, 115] aiming for a sensitivity of 
~ 10~ 10 pb is presented as well. One can see that Stage-2 detectors, like SuperCDMS, and 
eventually Stage-3 detectors, like XENON-1 ton, will be able to observe the signal from 
neutralino scattering off the nuclei in the entire HB/FP region and measure the respective 



6 Recent experimental results suggest a somewhat different value of the E term than the canonical value 
we assumed in this work. This affects the contribution from s-quark diagrams and can change our predictions 
for a si (Zip) by about a factor of three [ 108 1 . 
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Figure 14: Allowed region for spin-independent neutralino scattering cross section on proton 
plotted versus ■ The color code is the same as in Fig. [ll]. We also show the reach and projected 
reach of XENON-10, CDMS-II, SuperCDMS and XENON-1 ton detectors by solid, dashed, long- 
dashed and dotted lines, respectively. 



neutralino mass exhibiting the prominent complementarity to LHC. Observation of the 
signal in CDMS2 or a Stage-2 detector will clearly be an indication of the HB/FP region. 
These results can be used in conjunction, for example, with the LHC results for nig and rrih 
to determine Aq from Fig. 12, since the ambiguity of the HB/FP and the coannihilation 
region curves could be resolved using DD results. Also, positive result from DM searches 
will provide us with the neutralino mass, which, in turn, can be employed to extract Aq 
from the correlations of Fig. [D]. If the sensitivity of DD experiments would be further 
increased to ~ 10 -10 — 1CP 11 pb level, then individual curves from the stau-coannihilation 
pattern in Fig. 14 could be resolved and one could probe the Aq parameter in this region. 



4. Conclusion 

We have performed an updated scan of the CMSSM parameter space, taking into account 
the revised (lower) value of top quark mass, updated SUSY constraints from collider and 
low energy physics as well as crucial constraint on dark matter abundance from WMAP3. 
The scan was performed in the CMSSM parameter space for \i > 0, and for a plausible 
range of values for mo, mi/ 2 and |Ao| and for tan/? =5, 10, 50 and 53. 

We have demonstrated that taking account of the SUSY constraints, especially the 
dark matter abundance, strong correlations occur between the sparticle and Higgs masses. 
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The correlations between the light CP-even Higgs boson mass and SUSY particles could 
potentially allow determination of the sparticle spectra with a few percent accuracy, as- 
suming a theoretical control of rrih at the percent level. All correlations we have found are 
represented by separated narrow bands, two-fold solutions which exhibit focus point and 
h-funnel regions, or co-annihilation and A-funnel regions respectively. 

The correlations found among the sparticle and Higgs masses would also allow one 
to delineate the value of the trilinear coupling which can potentially be large. The large 
value of the trilinear coupling Aq is especially motivated by the present light Higgs mass 
constraints and the latest top-quark mass measurements. Since the top quark world- 
averaged mass went down to 170.9 GeV, the additional contributions from At and A), are 
even more important to increase rrih and thereby overcome the LEP2 limit. We have found 
that even for a parameter space with rriQ and \Aq\ as large as 5 TeV and 2 TeV respectively 
at Mgut, the lightest Higgs boson mass is limited by 122-123 GeV from above. 

We have found that some correlations, like between (tjia, rrih), ("M, rrig) and (rriA, rn^ ) 
could allow one to reveal the value of Aq and help to distinguish the DM motivated band 
of the parameter space mentioned above. Moreover, we have demonstrated striking com- 
plementarity between the LHC and direct dark matter detection experiments in resolving 
the SUSY mass spectrum and determination of SUSY parameters. Stage 2 experiments 
like SuperCDMS will be able to cover the entire HB/FP region and resolve the ambiguity 
of HB/FP and non-HB/FP bands for several important correlations. 

The correlations we have presented in this paper can be a useful input for SUSY global 
analysis fit as well as for experimental delineation of the SUSY parameter space. We would 
also like to emphasize that including the experimental input both from collider physics 
and from dark matter detection experiments would allow one to significantly improve the 
understanding of the SUSY spectrum and the underlying parameter space. 
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